Geometric magic numbers of sodium clusters: Interpretation of the melting behaviour 
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Putative global minima of sodium clusters with up to 380 atoms have been located for two 
model interatomic potentials. Structures based upon the Mackay icosahedra predominate for both 
potentials, and the magic numbers for the Murrell-Mottram model show excellent agreement with 
the sizes at which maxima in the latent heat and entropy change at melting have been found in 
experiment. 

PACS numbers: 61.46.+w,36.40.Mr 
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The melting of sodium clusters has been the subject 
of numerous recent studies. Much of this interest has 
arisen due to the availability of high-quality experimen- 
tal data, which has allowed detailed comparisons between 
theory and experiment. In particular, the Haberland 
group P, 0, Q have measured the caloric curves of mass- 
selected positively-charged sodium clusters with up to 
360 atoms, from which the melting temperature (T me jt) 
and latent heat can be extracted. The values of T me zt for 
these clusters are on average one third lower than that 
for bulk, and show variations of up to ± 50 K depending 
on the cluster size. There have been a number of theo- 
retical studies that, using different levels of theory, have 
investigated the origins of the size-dependence of T me i t 
0,0,0, 13 However, the peaks in the melting temper- 
ature do not seem to correlate either with the electronic 
or geometric shell closings of sodium clusters, and none 
of those theoretical studies have been able to provide a 
satisfactory explanation for the non-monotonic variation 
of T me i t . Significant progress was made in Haberland 
et a/.'s most recent paper, in which they observed that 
the energy and entropy changes on melting provide more 
structural insight into the system than T me u itself . In 
particular, these two quantities exhibit pronounced max- 
ima at certain 'magic numbers', some of which have a 
clear interpretation in terms of geometric structures, such 
as the Mackay icosahedra, whilst others remain unas- 
signed. Therefore, a systematic investigation of the ge- 
ometric structure of sodium clusters in this size range 
would be of great help in the identification of the struc- 
tures underlying these magic numbers. 

Previous work on the structure of sodium clusters has 
for the most part concentrated on clusters with less than 
60 atoms [Ifl[HHH[llH[li. By contrast, in 
this Letter we have attempted to locate the lowest-energy 
structures of sodium clusters for all sizes up to 7V=380 
using the basin-hopping global optimization method 
Such large sizes necessitate the use of a model potential, 
and we have considered two different forms for the in- 
teratomic interactions, namely the Gupta 

Hi 113 

and 

Murrell-Mottram (MM) [13. Ell E3| potentials. The MM 
potential has more parameters, has been fitted to a wider 
range of properties, and exhibits good transferability [22^ . 



Therefore, it is expected to be the more reliable of the 
two potentials, but it is also significantly more expensive 
to compute. The advantage of considering two potentials 
is that we can have greater confidence in those structural 
features that are common to both potentials. 

In Figure ^ we have plotted the energies of the puta- 
tive global minima for the two potentials, and Figure El 
shows the structures of some of the magic number clus- 
ters. The energies and coordinates for all the structures 
are available at the Cambridge Cluster Database [23| . 
For N < 57 the Gupta global minima have been previ- 
ously reported by Lai et al. jT3 |. 

The Haberland group found that for N < 100 many 
sodium clusters do not show a clear melting transition, 
but pass from solid to liquid without a pronounced la- 
tent heat Q. Na55 stands in contrast to this trend hav- 
ing a particularly high melting temperature, but Na7o 
and Nag2 also represent exceptions |9j. Both potentials 
exhibit a pronounced magic number at iV=55, which, 
as expected, corresponds to a complete Mackay icosa- 
hedron. Typically, there are two types of overlayer for 
growth on the surface of a Mackay icosahedron. The first, 
the Mackay overlayer, continues the face-centred-cubic 
(fee) packing of the twenty fee tetrahedra making up the 
Mackay icosahedron, and leads to the next Mackay icosa- 
hedron. By contrast, the second, the anti-Mackay over- 
layer, adds atoms in sites that are hexagonal close-packed 
with respect to the underlying fee tetrahedra. Typically, 
growth starts off in the anti-Mackay overlayer because of 
a greater number of nearest-neighbour interactions, but 
then switches to the Mackay overlayer because it involves 
less strain [ilEsT]. 

Interestingly, structures that do not adopt either of 
these over layers are prevalent for both potentials. The 
magic number at Nan, a possible explanation for the 
experimental feature at iV=70, provides a good example. 
Both potentials have the same C5 global minimum, where 
the five faces around the vertex of the 55-atom Mackay 
icosahedron are covered by a Mackay-like cap, but where 
both the overlayer and core have been twisted with re- 
spect to the ideal Mackay sites. This twist increases the 
coordination number of some of the surface atoms at the 
expense of increased strain and creates a structure where, 
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FIG. 1: Energy of the global minima found for the Gupta (upper panel) and MM (lower panel) potentials as a function of 
size. Energies are given relative to E ma gi C , which is a function fitted to the energies of the first four stronger magic numbers. 
0.2546iV 1/3 + 1.2134AT 2 / 3 - 1.15687V; E%M { c = -0.4788 + 0.5261iV 1/3 + 0.9852iV 2/3 - 1.1110JV. 



unlike both the anti-Mackay and Mackay overlayers, the 
surface consists entirely of {lll}-like faces. A similar 
structure is a magic number at N=92 and involves the 
covering of ten faces with a Mackay overlayer, which then 
undergoes a twist distortion, giving rise to a structure 
with T point group symmetry, instead of C^ v for the ideal 
Mackay geometry. These structures look like a hybrid of 
the 55-atom and 147-atom Mackay icosahedra, because 
they have triangular {111} faces of sizes corresponding 
to both the smaller and larger Mackay icosahedra. 

The Gupta potential does not exclusively exhibit struc- 
tures based on these twisted icosahedra in this size range. 
For example, at N=8l a structure where eight faces are 
covered by an anti-Mackay overlayer is the global min- 
imum. This difference between the two potentials be- 
comes more prominent at larger sizes. For example, there 
is a feature at N » 116 in the experimental results that 
has been interpreted in terms of a Mackay structure with 
15 of the faces of the underlying icosahedron covered @. 
Both potentials have clear features near to this size. For 
the MM potential, there is a magic number at iV=116 and 
at this size there are two minima with almost the same 
energy. The second-lowest minimum corresponds to a 
twisted form of the structure suggested by Haberland et 
al., and the lowest-energy isomer is based on a 116-atom 
Ino decahedron but with the central ring of atoms twisted 



to remove any {100} faces. By contrast, the Gupta po- 
tential has a magic number at iV=115 that corresponds 
to an Ih structure with a complete anti-Mackay overlayer. 
This is an unusual feature, since the anti-Mackay over- 
layer is usually observed during the initial growth on an 
icosahedron |2J], but not when that overlayer is nearly 
complete. Moreover, this structure is very high in energy 
for the MM potential. 

Experimentally, Nai47 is a prominent magic number, 
and, again as expected, the Mackay icosahedron is the 
global minimum at this size for both potentials. However, 
for the Gupta potential a more stable structure can be 
obtained by removing the twelve vertex atoms, giving rise 
to a magic number at 7V=135 (Fig. [I}. This feature is in 
clear contradiction with experiment. 

For growth on the 147-atom Mackay icosahedron, the 
differences between the results for the two potentials be- 
come even greater. For the MM potential structures 
based upon the twisted icosahedra continue to predom- 
inate. However, the Gupta potential initially exhibits 
structures with a Mackay overlayer, and then switches 
over to an anti-Mackay overlayer near to the completion 
of this overlayer at 7V=267 . 

The MM potential exhibits prominent magic numbers 
at 7V=178, 216, 232 and 258, with weaker features at 
184, 190, 200, 206, 222 and 238. These structures cor- 
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FIG. 2: A selection of sodium clusters that show enhanced stabilities for the (a) Gupta and (b) Murrell-Mottram potentials. 
Each structure is labelled by the number of atoms and its point group symmetry. 



respond to covering successive faces of the 147-Mackay 
icosahedron with Mackay-like overlayers, but where the 
core and surface again undergo a twist distortion. The 
178-, 216-, 232- and 258-atom structures are equivalent 
to the 71-, 92-, 101- and 116-atom structures described 
above and correspond to covering all the faces surround- 
ing 1, 3, 4 and 6 vertices of the underlying icosahedron. 
These features are in good agreement with the experi- 
mental results, which have clear features at iV=178 and 
216, and a smaller sub-peak at iV=184. No experimen- 
tal features have yet been identified at A^=232 and 258. 
However, in this size range the data is sparse, and the 
error bars are of similar magnitude to the size variation 
of the properties. Therefore, it would be very interesting 
if further experiments were conducted at these sizes to 
examine the predictions of the MM model. 

Interestingly, Haberland et al. suggested undistorted 
Mackay structures to explain the magic numbers at 
N=178 and 216 9J. However, it is more usually found 
that more stable structures are possible, when the five- 
coordinate atoms at the corners of the added triangular 
faces are not occupied. For example, this leads to magic 



numbers at iV=173 and 213 for Lennard-Jones clusters 
[i^l . The twist distortion of the icosahedra provides a 
possible explanation for this difference in magic numbers. 
As a consequence of the distortion, the coordination num- 
ber for the corner atoms increases from five to six, making 
it more favourable for these sites to be occupied. 

The magic numbers for the Gupta potential are com- 
pletely different in this size range, because of the pref- 
erence for both undistorted icosahedral structures and 
empty vertex sites. The magic numbers at A r =166, 186, 
201, 216 and 241 are all due to structures with a Mackay 
overlayer. If it were favourable for the six-coordinate ver- 
tices to be occupied, these magic sizes would instead be 
at A/=173, 196, 213, 230 and 258. Only if five-coordinate 
sites were also occupied would the magic numbers be 178, 
200, 216, 232 and 258. Analogous to the particular sta- 
bility of Nans in the growth of the third shell, there is an- 
other magic number at N—267 whose structure involves 
a complete anti-Mackay layer without vertices. Closeby 
(N = 268), Haberland et al. found a well-structured pho- 
toelectron spectrum, but they attributed this feature to 
the existence of an electronic shell closing rather than to 
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high point group symmetry |9j. Furthermore, this com- 
plete anti-Mackay icosahedra again lies very high in en- 
ergy for the MM potential. 

As for the third shell, the complete Mackay icosahe- 
dron is not a magic number for the Gupta potential, 
but instead an icosahedron with twelve missing vertices 
is more stable, displacing the magic number to iV=297. 
The MM potential still predicts the magic number to be 
at 7V=309, but the difference in stability between the 
297- and 309-atom structures is much smaller. Indeed, 
at N w 360 structures with missing vertices actually be- 
come more stable. The similar behaviour of the two po- 
tentials suggests that the loss of vertex atoms is a robust 
structural feature for sodium; the two potentials only dif- 
fer in the size at which this effect first appears. These 
results suggest that a plausible explanation of the absence 
of an experimental magic number at iV=309, but the ap- 
pearance of a feature at iV=298, is the greater stability of 
a Mackay icosahedron that has lost its vertices. However, 
Haberland et al. found that the measured photoelectron 
spectrum for Na298 is not compatible with such a struc- 
ture Furthermore, on measuring the photoelectron 
spectrum of Na3og as a function of temperature, Haber- 
land et al. found that a structural transition occurred at 
about 40 K below melting |9j . Parallel tempering simu- 
lations using the Gupta potential, however, did not show 
evidence of any transitions prior to melting for Na3og. 
Therefore, our results are unable to offer a structural ex- 
planation compatible with all the experimental findings 
associated with the completion of the fourth icosahedral 
shell, and their origin remains somewhat mysterious. 

Finally, for growth of the fifth icosahedral shell, the 
same patterns continue, i.e. Mackay over layers for the 
Gupta potential and twisted icosahedra for the MM po- 
tential. In this size range, experiments predict a peak 
at 7V=360, for which the MM magic number at 7V=357 
offers a possible explanation. 

In summary, our results support the conclusions de- 
rived from Haberland et al.'s recent analysis of the melt- 
ing behaviour of sodium clusters — namely that the clus- 
ters in this size range are predominantly icosahedral. Of 
the two potentials we have considered, the structures ob- 
tained for the MM potential appear more reliable and 
the magic numbers are in excellent agreement with ex- 
periment. In particular, we suggest that the experimen- 
tal features at sizes intermediate between the complete 
Mackay icosahedra are due to icosahedral structures with 
a Mackay overlayer, but where both the core and over- 
layer undergo a twist distortion to give structures that 
have only {lll}-like faces. It is noteworthy that such 
features cannot be captured by a pairwise-additive inter- 
atomic potential. 
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